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ABSTRACT 

' Very high-energy y-ray emission from PKS 0447-439 was detected with the H.E.S.S. Cherenkov telescope array in December 2009. This blazar 
» I , is one of the brightest extragalactic objects in the Fermi Bright Source List and has a hard spectrum in the MeV to GeV range. In the TeV range, a 
& . photon index of 3.89 ± 0.37 (stat) ± 0.22 (sys) and a flux normalisation at 1 TeV, OiTeV = (3.5 ± l.l(stat) ± 0.9(sys)) x lO^cirrV'TeV -1 were 
found. The detection with H.E.S.S. triggered observations in the X-ray band with the Swift and RXTE telescopes. Simultaneous UV and optical data 
from Swift UVOT and data from the optical telescopes ATOM and ROTSE are also available. The spectrum and light curve measured with H.E.S.S. 
are presented and compared to the multi-wavelength data at lower energies. A rapid flare is seen in the Swift XRT and RXTE data, together with a 
flux variation in the UV band, at a time scale of the order of one day. A firm upper limit of z < 0.59 on the redshift of PKS 0447-439 is derived 
from the combined Fermi-LAT and H.E.S.S. data, given the assumptions that there is no upturn in the intrinsic spectrum above the Fermi-hAT 
energy range and that absorption on the Extragalactic Background Light (EBL) is not weaker than the lower limit provided by current models. The 
spectral energy distribution is well described by a simple one-zone Synchrotron Self-Compton (SSC) scenario, if the redshift of the source is less 
than z< 0.4. 
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1. Introduction 

With only few exceptions, the Active Galactic Nuclei (AGN) so 
far detected in very high-energy y-ray s (VHE; E > 100 GeV) 
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belong to the class of blazars and are thought to have their jets 
closely aligned to the line of sight, which leads to an amplifi- 
cation of their apparent luminosity through relativistic beaming. 
While flat-spectrum radio quasars (FSRQ) account for an impor- 
tant fraction of sources detected at ene rgies in the MeV to GeV 
range (>40% of blazars of known tvpe. lAckermann et al.ll201 iK 
the population of blazars detected at VHE is largely dominated 
by BL Lac-type objects — especially hi gh-frequency peaked B L 
Lac objects (HBL), as defined by Padovani & Giommi ( 1995). 

A general difficulty in the study of BL Lac objects is that 
their redshifts are often unknown or not well determined, since 
spectral features from the host galaxy become undetectable due 
to the dominance of the non-thermal emission from the jet. 
Therefore, about 50% of the BL Lac objects detected with 
Fermi-LAT (lAckermann et al.ll201 ll) have an unknown redshift, 
and this is also the case for PKS 0447-439. The so urce was 
first discovered in the radio band ([Large et al.lll98l|) and has 
since been detected in several observation s in the radio, in- 
frared, optical , ultraviolet and X-ray bands (IWhite et al. 1994t 
Gregory et al.1 1 19941 lLamptqn et al.1 119971: Craig & Fruscione 
1997t lHaakons en & Rutledge 2009). It was identified as a bright 
B L Lac object b y Perlm aTTet al.l (1 1 998 ) and classified as an HBL 
bv lLandt & Bignall d2008l) . based on the ratio of the radio core 
luminosity at 1 .4 GHz over the X-ray luminosity at 1 keV. 

A first evaluation of its redshift of z= 0.107 
dCraig & Fruscionel 19971) was based on a mis-identification with 
a Seyfert 1 galaxjQ- Subsequent observatio ns with the CTIO 4m 
telescope led to an estimation of z=0.205 (Perl man et al.l fl998). 
However, this claim was based only on a very weak spectral 
feature that the authors identifed as the Ca II line in an otherwise 
featureless spectrum. This result was not confirmed by a more 
recent atte mpt at determining the redshift of the source. The 
analysis by lLandt & Bignalll (120081) yielded only a lower limit 
of z>0.176, based on t he photometric method described by 
Piranomo nte et al.l (120071) . which was applied to a featureless 
spectrum resulting from observations with the CTIO and NTT 
telescopes. 

A recent analysis of observations of the source, which were 
carried out with the CTIO and NTT i n 2007, puts a very high 
lower limit of z > 1.246 on its redshift (Landt 2012). This result 
is based on weak absorption lines interpreted as the Mg II /12800 
doublet. Such a high redshift is very difficult to reconcile with 
the detection of y-rays from the source with energies of several 
TeV, given the absorption by the extragalactic background light 
(EBL). A more likely explanation is that the feature in the opti- 
cal spectr a corresponds t o atmospheric absorption. The most re- 
cent data dPita et alJl2012l) . taken with the X-Shooter telescope 
(ESO/VL T), do not confirm the high redshift value suggested 
bv lLandtl d2012l) . but show instead a featureless spectrum apart 
from atmospheric absorption lines, at wavelengths that coi ncide 
exactly with the putative Mg II doublet suggested by lLandtl 
d2012l) . This has also been confirmed subsequ ently by obser- 
vation s with the Clay Magellan II telescope ( IFumagalli et al. I 
l20T2h . 

While the optical spectrum of PKS 0447-439 is featureless 
and seems strongly dominated by non-thermal emission from the 
jet, observations with AT CA in the radio band sh ow an extended, 
lobe-dominated source dLandt & Bignalll [2008). The emission 
from the jet, although strong at optical energies, is much weaker 
at radio energies, possibly due to synchrotron self-absorption. 

High energy (HE) emission in the MeV to G eV range from 
the source has been detected with Fermi-LAT dAtwood et al. I 



private communication with H. Landt and M. Veron-Cetty 



2009). PKS 0447-439 is associated with 0FGL J0449.7-4348, 
one of the bright est blazars in the Fermi Bright Source List 
dAbdo et alj 120091) . and has a hard spectrum in the HE band, 
which has made it a prime target for observations in the VHE 
range. In 2009, VHE emission from the source was discov- 
ered with the H.E.S.S. (H igh Energy Stere oscopic System) 
Cherenkov telescope a rray dRaue et al. 120091) and preliminary 
results were presented (IZech et al.ll201 ll) . 

In the present work, the final results of the H.E.S.S. analysis, 
using an updated analysis chain, are presented, together with si- 
multaneous data at other wavelengths. Constraints on the emis- 
sion processes and on the redshift of the source are discussed. 
Details on the VHE data analysis are provided in Section|2] The 
discovery of VHE emission triggered observations of the source 
with the RXTE and Swift space telescopes in the X-ray, UV and 
optical bands. All the available multi-wavelength (MWL) data, 
including also optical data from the ATOM and ROTSE tele- 
scopes, situated on the H.E.S.S. site, are presented in Section[3] 
The spectral energy distribution (SED) and the temporal flux 
evolution in different energy bands are discussed in Section [4] 
The comparison of the spectra measured with ferm/-LAT and 
H.E.S.S. is used to derive a firm upper limit on the redshift of the 
source. An interpretation of the SED following a standard syn- 
chrotron self-Compton (SSC) scenario is suggested in Section|5] 
for different assumptions on the redshift of the source. 

Throughout the paper, a concordance cosmology model with 
Hq = 70km s Mpc -1 , Om = 0.3, and £2a = 0.7 is assumed. 
Unless otherwise indicated, all errors are statistical and are given 
at the lcr (standard deviation) confidence level. 



2. H.E.S.S. data analysis 

A total of 13.5 hours of good quality data were taken between 
November 2009 and January 2010 at the source position of 
PKS 0447-439 with the Cherenkov telescopes of the H.E.S.S. 
array, located in the Khomas highland in Namibia. Data taken 
when less than three of the four telescopes were fully operational 
have been rejected. Most of the data presented here are from 
December 2009. Only 0.9 hr of observations from November 
2009 and 0.4 hr from January 201 were kept in the final data set 
after standard quality selection (Aharonia n et al . 2006). Zenith 
angles of the observations range from 20° to 26°. 

Data were analysed using the model analysis 
dde Naurois & Rollandll2009l) . which yields a more than twice 
as high si gnificance for this s ource than a standard Hillas-type 
analysis (Ahar onian et al.l2006l) . mainly due to a better rejection 
of the cosmic-ray background. For the given range of zenith 
angles and the minimum requirement of 60 photoelectrons 
per image used in the analysis, the analysis yields an energy 
threshold of about 220 GeV. All results were cross-checked 
with independent ana lysis procedures and calibration chains 
dAharonian et al.ll2006l) . The data yield a strong VHE signal at a 
statistical significance of 15.2<r. The flux measured above 250 
GeV was ~3% of the flux from the Crab Nebula measured with 
H.E.S.S above the same threshold. The best fit position of the 
VHE y-ray excess (072000 = 4 h 49 m 28.2 s ± 1.4 s (stat) + l s (sys), 
5/2000 = -43°50'12" + 15" (stat) +18" (sys)) is com patible 
with the nominal position dVeron-Cettv & Veron ll2010T) of the 
source (072000 = 4 h 49 ra 24.7 s , c5 y2 ooo = -43°50'9"). 

The distribution of the y-ray excess as a function of the 
squared angular distance from the nominal source position 
(Fig. Q} has been determin ed using the "reflected background" 
method (Berg e~et"a"fl l2007|) . From the on- and off-region, 336 
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and 1439 events were detected, respectively, with a relative nor- 
malisation factor between the regions of 0.083. The distribution 
of the on-source events is consistent with that expected from a 
point source. 




Fig. 1. Distribution of the squared angular distance of y-ray 
events from the nominal source position. The points represent 
the distribution for on-source events. The normalised off-source 
event distribution is described by the solid histogram. The verti- 
cal dashed line indicates the on-source integration region of 0.01 
square degrees used in the analysis. 



An average power-law spectrum measured over the whole 
observation period has been extracted using forward fo lding with 
a maximum likelihood optimisation dPiron e t al. 2001). The de- 
rived VHE spectrum is very soft, with a photon index T = 3.89 
+ 0.37 (stat) ± 0.22 (sys) and a flux normalisation at 1 TeV of 
(Dixev = (3.5 ± 1.1) x lO^cirrV^eV- 1 . The equivalent X 2 of 
the fit is 36 for 26 degrees of freedom (d.o.f.) with a chance 
probability of 0.1. The systematic uncertainty in the flux nor- 
malisation is about 25%. There is no significant indication for a 
spectral break or curvature in the observed spectrum. 

The nightly binned integrated flux above 250 GeV has been 
extracted using the measured photon index (cf. Table Q] and 
Fig.|2]i. No significant signal is seen during the short observation 
time in November 2009 and January 2010. A fit of a constant to 
the nightly binned fluxes detected only during December 2009, 
weighted by the statistical uncertainties, yields an average flux 
level of (4.7+1.0)xl0- 12 cirT 2 s _1 with a X 2 of 16.3 for 9 d.o.f. 
and a chance probability of 0.06. 



When including the data from November and January in the 
fit (in nightly bins), a lower average flux of (3.3 ± 0.8) x 10~ 12 
cm -2 s _1 with a x 1 of 22.2 for 11 d.o.f. (chance probability 
of 0.02) is found. T he fractional RMS variability amplitude 
dVaughan et al . 2003) over the complete observation period is 
F V ar = 0.82±0.20. This is marginal evidence that the VHE flux 
was higher during December 2009 compared to the previous and 
the following month, although the statistics for these months are 
very limited due to the short duration of good-quality observa- 
tions. 



<> 25 
E 
o 

? 20 
o 

?15 
a> 

a 

o 10 



~ 5 



3E4 



55140 55150 55160 55170 55180 55190 55200 55210 



MJD 



Fig. 2. H.E.S.S. light curve from November 2009 to January 
2010 in nightly bins with a constant fitted to the December data 
(dashed line) and to the whole dataset (dash-dotted line). 
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Table 1. Nightly averaged fluxes measured with the H.E.S.S. 
array. For each night, the start to and end ti of observations (in 
MJD), the live time At (in seconds) and the integrated y-ray flux 
above a threshold of 250 GeV (in units of 10~ 12 cm~ 2 s~') are 
shown. 



3. Multiwavelength data 

3.1. High Energy data from Fermi-LAT 

PKS 0447-439 has been detected in all three Fermi-LAT point- 
source catalog ues and is associated with the source 2FGL 
J0449.4-4350 dAckermann et al.ll201 ll) in the latest catalogue. A 
fit of a constant to the light curve from the Fermi-LAT two-year 
point-source catalogue yields an average photon flux between 
0.1 GeV and 100 GeV of (7.75 ± 0.24) x 10~ 8 cirT 2 s _1 . The 
energy flux in this range is quoted as (1.36 + 0.06) x 10~ 10 erg 
cm -2 s _1 . The photon index, derived from the likelihood analysis 
for 0.1 GeV to 100 GeV, has a value of 1.86 ± 0.02, indicating a 
hard spectrum. The source is seen to be variable (variability in- 
dex of 91.9), with the maximum peak occurring on MJD 54757 
(October 10, 2008). The monthly averaged flux around this peak 
between 0.1 GeV and 100 GeV is (1.56 + 0.16) x 10~ 7 cirT 2 s _1 . 

Fermi-LAT data have been analysed between MJD 54682 
(August 4, 2008) and MJD 55819 (September 15, 2011), in an 
energy range from 0. 1 GeV to 200 GeV. The Fermi-LAT analysis 
was performed on the publicly available data, using the binned 
likelihood method (cf . lAckermann et alj[201 ll) from the Science 
Tools package, version v9r23plQ provided by the Fermi-LAT 
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Fig. 3. H.E.S.S. spectrum of the source extracted from the com- 
plete data set. The differential flux points (with an excess signif- 
icance of lcr) and the 68% confidence band for the best fit are 
shown in the upper panel, as a function of true energy. The arrow 
represents an upper limit at the 99% confidence level. The resid- 
uals of the maximum likelihood fit are given in the lower panel 
as a function of the reconstructed energy, together with a fit of a 
constant (dashed line). 



collaboration. Source class events are considered in a circular 
region of interest with a radius of 10° around the nominal posi- 
tion of PKS 0447-439, using the P7V6.SOURCE instrumental re- 
sponse functions. A cut was applied on the zenith angle with re- 
spect to the Earth (<100°), and on the rocking angle (<52°). All 
the objects included in the 2FGL within 15° were included in 
the model reconstruction of PKS 0447-439. The isotropic model 
iso_p7v6source was used to account for both the extragalac- 
tic diffuse emission and residual instrumental background, while 
the spatial template gal_2yearp7v6_v0 was used to account for 
the contribution from the Galactic diffuse emission. 

This analysis yields a photon index of 1.86 ± 0.02 and an 
integrated energy flux between 0.1 GeV and 200 GeV of (1.45 
+ 0.05) x 10~ 10 erg cirT 2 s _1 . A long-term light curve with a 
binning of 90 days has been extracted, assuming a power-law 
spectrum with a fixed index of 1.86 (Fig. |4j. A second analysis 
has been carried out with a maximum energy of 100 GeV, as 
used for the 2FGL, yielding an integrated energy flux between 
0.1 GeV and 100 GeV of (1.26 ± 0.05) x 10" 10 erg cirT 2 s _1 
and a photon index of 1.85 + 0.02, in agreement with the results 
from the 2FGL. The estimated systematic uncertainty on the flux 
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Fig. 4. Fermi-LAT long-term lightcurve of PKS 0447-439 
showing the photon flux integrated between 0. 1 GeV and 200 
GeV in intervals of 90 days. The horizontal error bars indicate 
the duration of the intervals. The period chosen for the extraction 
of the Fermi-LAT spectrum is indicated by the large blue band. 
The thin, striped band corresponds to the H.E.S.S. observational 
period in December 2009. 



in these analyses is 10% at 0.1 GeV, 5% at 0.5 GeV and 10% at 
10 GeV and above (cf. lAckermann et al.ll201 lb . 

In the Fermi-LAT long-term light curve (cf. Fig.|4}, a period 
of relatively constant flux was identified between August 7, 2009 
(MJD 55050) and December 20, 2009 (MJD 55185), an interval 
which includes the nights of H.E.S.S. observations of the source 
in 2009. For this selected period, an averaged spectrum with a 
power-law shape was extracted with a photon index of 1.85 + 
0.05 (included in Fig. |7J. Using the gtlike tool and assuming a 
power-law shape for the source spectrum, the test statistic (TS) 
of the likelihood analysis of PKS 0447-439 is 1 143, correspond- 
ing approximately to a 34<x detection. The TS is 1143.43 for 
the power-law fit, 1143.66 for a log-parabolic fit, and 1146.23 
for a fit with a broken power law. Comparing the log-parabolic 
model with the power law yields a log likelihood ratio (LLR) of 
0.30. For the comparison between the broken power-law and the 
power-law model, the LLR is 2.96. Thus, spectral shapes more 
complex than a power law do not result in an improvement of 
the likelihood for PKS 0447-439. 

For the period of constant flux (MJD 55050-55185), an in- 
tegrated flux between 0.1 GeV and 200 GeV of (1.08 ± 0.09) 
x 10 -7 cirT 2 s _1 was found. A fit of a constant to the selected 
period of the Fermi-LAT light curve (in bins of 3 day) yields 
a chance probability of 0.98 Of 2 /d.o.f. = 27.3 / 44). The most 
energetic photon detected within the 95% point-spread function 
(PSF) containment radius was found at 4.6 arcminutes from the 
nominal source position, with an energy of 95.5 GeV. For the full 
period (MJD 54682-55819), the most energetic photon within 
the PSF of the source has an energy of 147.5 GeV, while the en- 
ergy resolution at this energy is about 12%. A part of the Fermi- 
LAT light curve, binned in intervals of three days, is included in 
Fig.S 

3.2. X-ray data from RXTE and Swift XRT 

The discovery of a VHE signal from PKS 0447-439 triggered 
observations with the RXTE and Swift space telescopes. In 
December 2009, 11 pointings were taken with the PCA units 
of RXTE dJahoda et alJ Il996l) . The STAND ARD2 data were 
extracted using the HEASOFT 6 . 9 analysis software package 
provided by NASA/GSFC, and filtered using the RXTE Guest 
Observer Facility (GOF) recommended criteria. Only signals 
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from the top layer of the PCA were selected. The faint back- 
ground model was applied. Due to low statistics, the data from 
several consecutive pointings had to be added together to allow 
the extraction of energy spectra (cf. Table [2}. The spectra were 
extracted between 3 and 12 keV using XSPEC v . 1 2 . 6 . 8, with 
a fixe d column density of Nu = 1 .24 x 10 20 cirT 2 dKalberla et alj 
I2005h . 

A power-law model was used for spectral fitting. No signif- 
icant improvement was noted when using a broken power law. 
The parameters of the spectra for the different pointings are 
listed in Table [2] and the light curve is included in Fig. [5] and 
will be discussed in Section l4~Tl 

In the soft X-ray band, the source was observed with Swift 
XRT dBurrows et al.l 120051) in 7 pointings taken in December 
2009. The analysis of the XRT data (in photon-counting mode) 
was carried out with the HEASOFT 6 . 9 standard tools. Source 
events were extracted within a circle with a radius of 20 pixels 
(~47 arcseconds), while background events were extracted from 
both an annular region around the source (between 30 and 80 
pixels) and a nearby source-free circle with a radius of 80 pixels. 
It was found that in some pointings pile-up was present. In these 
cases, the internal part of the PSF was excluded in the event 
selection and data were extracted in an annular region with an 
outer diameter of 30 pixels, while the inner radius was chosen 
by fitting the radial distribution with a King profile. Using the 
grppha tool, bins were grouped to yield at least 30 counts per 
bin, and the data were fitted with a single or broken power law, 
with the same fixed column density as given above. Only points 
in the intervals from 0.3 to 0.45 keV and from 0.6 to up to 6 keV 
were included in the fit to suppress known systematic effects at 
intermediate energies (Campana et al. 2006). 

The results from the power-law and broken power-law fits 
were compared using an F-test, where a broken power law with 
a break at around 1 keV was preferred at a probability of >99% 
for two pointings, as indicated in Table [3] The hypothesis of ad- 
ditional, intrinsic absorption was also tested and was not found 
preferable over the power laws with Galactic absorption. The 
light curve of the soft X-ray flux is included in Fig.[5]and will be 
discussed in Section l4~T1 

The data taken with Swift XRT and RXTE PCA are not si- 
multaneous, with the exception of MJD 55184 (Dec. 19, 2009), 
where an overlap of about 0.7 ks occurs. For this night, data from 
the first pointing with PCA have been jointly analysed with data 
from the second pointing with XRT. A fit of the combined spec- 
tra with a power law yields a photon index of 2.31 ± 0.03 and a 
deabsorbed energy flux of (1.16+q jj) x 10~ n erg cirT 2 s _1 be- 
tween 2 and 10 keV for a fixed N H = 1.24 x 10 20 cirT 2 . This fit 
has a reduced^ 2 of 0.91 for 72 d.o.f. A fit with a broken power 
law does not yield a significant improvement. 

3.3. Ultra-violet and optical data from Swift UVOT, ATOM and 
ROTSE 

Simultaneously to the Swift XRT pointings, observations in the 
ultra-violet (UV) a nd optical bands we re made with the Swift 
UVOT instrument (Ro ming et al.ll2.Q05h . Data were taken with 
the six available filters (v, b, u, UVW1, UVM2, UVW2 - in order 
of increasing frequency). UVOT counts were extracted within 
a radius of 5 arcseconds centered on the source. Flux densities 
have been extracted with uvotmaghist vl.l and have been mul- 
tiplied with the central wavelength of each filter dPoole et al.1 
20081). Galact ic extinction has been corrected for, following 
Roming et alJ (120091) . with an extinction coefficient Eb-v = 



0.014 (Schleg el etal1 fT998). The light curves are included in 
Fig.|5]and nightly averaged fluxes are listed in Table|4] 

Apart from the data from the Swift UVOT instrument, 
optical data are also av ailable over longe r intervals from 
the ground-based AT OM (lHauser et al J 120041) and ROTSE-IIIc 
(lAkerlof et all 120031) telescopes, both located on the H.E.S.S. 
site. 

Data from the 75 cm ATOM telescope are available in the B 
and R bands from September 2009 to December 2012. During 
the H.E.S.S. campaign, only few observations were taken in the 
B band and some of them were of insufficient quality to yield 
acceptable flux measurements due to the presence of clouds. A 4 
arcsecond radius of aperture was used to extract the data. In the R 
band, the ATOM light curve (Pig. [5j shows a slow flux decrease 
in December 2009. In the 2009 ATOM data and up to mid 2010, 
the magnitude increased by about 0.7 mag. The source was still 
in a relatively high flux state in December 2009 as compared to 
the following months. On longer time scales, variations of about 
1 mag were observed, and the source reached again about the 
same flux level as in December 2009 by the end of 2010 and by 
the end of 2011. 

ROTSE has a wide field of view (1.85° xl.85°) and is oper- 
ated without filters. A relative R magnitude is derived by com- 
parison of_fhe_Jn^umental magnitude with the USNO cata- 
logue (lAkerlof et al.ll2000l) . The ROTSE light curve, included 
in Fig. [5] shows the same slow flux decrease during December 
2009 as seen in the ATOM data. Long-term optical data taken 
with ROTSE between 2005 and 2010 confirm that the optical 
flux was decreasing, but in a relatively high state in December 
2009. 

The data from both instruments have been corrected for 
Galactic absorption as described above. 

4. Results 

4.1. Flux evolution in different wavebands 

Figure [5] shows a comparison of the flux evolution in different 
energy bands during the time period where data from H.E.S.S. 
are available (November 2009 to January 2010). The first panel 
shows the H.E.S.S. light curve above 250 GeV in nightly in- 
tervals, with negative flux estimate s replaced by upper l i mits a t 
the 95% confidence level following Feldman & Cousins ( 1998). 
Details on the flux evolution in the H.E.S.S. data were given in 
Section [2] 

Over the same period, the Fermi-LAI light curve between 
0. 1 and 200 GeV, binned in intervals of three days, is shown in 
the second panel. (A light curve with one-day binning suffers 
from very large statistical uncertainties.) The fractional RMS 
variability amplitude F var = 0.21 ± 0.22 does not point to any 
significant intrinsic variability during this period. However, the 
fit of a constant to the light curve has a^ 2 of 75.3 for 24 d.o.f., 
corresponding to a chance probability of less than 10~ 6 , which 
indicates a deviation from a constant flux after the stable period 
between MJD 55050 and MJD 55185 (cf. Sec. ED, notably due 
to the low fluxes on MJD 55 195 and MJD 55207. When compar- 
ing the H.E.S.S. and Fermi-LAI light curves, it can be seen that 
the absence of a positive excess in the H.E.S.S. data in November 
2009 is reflected in the absence of a signal in the Fermi-LAI flux 
for the same night where the H.E.S.S. data were taken. 

Rapid variability on nightly time scales is clearly detected 
in the X-ray data (third panel). Fits to a constant of the Swift 
XRT and RXTE PCA light curves yield chance probabilities be- 
low 10~ 15 . The fractional RMS variability amplitude for the data 
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Table 2. Hard X-ray spectra measured with RXTE/PCA. The columns provide for each pointing the start to and end ti of observations 
(in MJD), the live time At (in seconds), the count rate (per second), the integrated deabsorbed flux between 2 and 10 keV (in units of 
10~ n erg cm 2 s _1 ), the photon index, reduced^ 2 and the number of degrees of freedom. Data from the observations in the second 
to fourth row and in the sixth to ninth row have been summed up to extract spectra with improved statistics. 
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Table 3. Soft X-ray spectra measured with Swift XRT. The columns provide for each night the start to and end ti of observations 
(in MJD), the live time 6t (in seconds), the count rate (per second), the integrated deabsorbed flux between 0.3 and 4 keV (in 10" 11 
erg crrT 2 s~'), the photon index below the break energy, the break energy (in keV), the photon index above the break energy, the 
reduced x 1 an d the number of degrees of freedom. For nights where no break energy is given, Fi is the photon index for a power 
law over the whole energy range. 
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Table 4. Nightly averaged fluxes in the optical and ultra-violet bands measured with SWIFT UVOT. The date of the nights is given 
in the first column (in MJD). The remaining columns provide the deabsorbed flux densities for the different filters (v, b, u, UVW1, 
UVM2, UVW2), multiplied with the central wavelengths (5468 A, 4392 A, 3465 A, 2600 A, 2246 A, 1928 A) in units of 10" 11 
erg cirT 2 s _1 . Errors represent the RMS spread between points for nights where more than one measurement was available (marked 
with an asterisk *) and the statistical uncertainty in the single measurement otherwise. 



from XRT is F var =0.22 ± 0.01 and for PCA F var =0.32 + 0.05. 
In the hard X-ray band (2- 1 keV), data from RXTE PCA show a 
flare on MJD 55 1 85 (December 20). The flux drops to about half 
its peak value on the following night. As far as can be concluded 
from the limited statistics available, the flare is not accompanied 
by any significant spectral variation. In the soft X-ray band (0.3- 
4 ke V), the Swift XRT light curve exhibits a flare on MJD 55183 
(December 18, 2009), where the flux reaches a level about twice 
as high as its lowest value in this observation period. No sig- 



nificant correlation is found between the photon indices and the 
integrated fluxes for the different pointings. The minimum vari- 
ability time scale in the X-ray bands ha s been estimated in terms 
of the flux doubling time T2, following Zha ng et al.l {,1999). The 
smallest flux doubling time was found in the XRT light curve 
with a value of T2 = 0.87 + 0.24 days (after discarding results 
with uncertainties above 30 %), indicating flux variations on the 
order of one day. 
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Fig. 5. Multi-wavelength light curves for the time period from November 2009 to January 2010. The first panel shows the integrated 
flux above 250 GeV seen with H.E.S.S., while the Fermi-LAT flux between 0.1 and 200 GeV is presented in the second panel. 
Horizontal error bars indicate the duration of the time intervals over which fluxes were averaged. Upper limits are given at the 95% 
confidence level (CL) in the H.E.S.S. and Fermi data. The third panel shows the X-ray fluxes between 2 and 10 keV from RXTE 
(open circles, y-axis on the left) and between 0.3 and 4 keV from Swift XRT (filled circles, y-axis on the right). UV data taken with 
Swift UVOT with the u, UVW1, UVW2 and UVM2 filters are shown in the fourth panel. The last panel includes optical data from 
the b and v filters from UVOT, in the B and R bands from ATOM and in the R band from ROTSE. For the UV and optical data from 
UVOT, flux densities have been multiplied by the central wavelengths of the corresponding filters. 



It is important to note that none of the pointings of the two 
X-ray telescopes were exactly simultaneous, with some overlap 
occurring only during one pointing on MJD 55184 (see Fig. |6). 
The flares seen in the two X-ray bands might be part of the 
same event, but this cannot be concluded with certainty. No cor- 
responding flare can be seen in the Fermi-LPtX data, but flux 
variations of the order seen in the X-ray band might well be hid- 
den in the statistical uncertainties. In the VHE band, no data are 
available during the night of the Swift XRT flare. 

The data taken with Swift UVOT in the UV band (fourth 
panel of Fig. [5]) show significant variation in the u- and UVW1- 
filters, when fit with a constant flux (chance probabilities < 
10 11 ; F var = 0.08 + 0.01 for the U-filter and F var = 0.09 + 0.01 
for the UVWl-filter). Variation in the UVW2- and UVM2-filters 
is only marginally significant. A variation is seen in the flux evo- 
lution that is simultaneous to the flare detected with Swift XRT. 
When looking at the nightly averaged UV fluxes, the flux in- 
crease in the UVOT data between MJD 55182 and MJD 55183 
is not significant, but there is evidence of a flux decrease of about 



10% between MJD 55183, where the X-ray flare is seen, and the 
following night (see Fig.|6]and Table[4]l. 

In the optical band (last panel of Fig.|5]l, flux variations dur- 
ing the period where the X-ray flares are seen are not statistically 
significant (with a 2% chance probability for a fit of a constant 
to the light curve in the UVOT v-filter), but there is an indication 
for a flux decrease between MJD 55183 and the following night 
in the v- and b-filters observed with Swift UVOT. This is not 
visible in the R band, covered with ATOM and ROTSE. ATOM 
data in the B band suffer from large uncertainties during several 
nights, since these observations are affected by clouds, and no 
significant variation was detected. 

4.2. Spectral energy distribution 

The spectral energy distribution of PKS 0447-439 is presented 
in Fig.|7](the models are discussed in Section|5]l. The VHE spec- 
trum is not corrected for absorption by the EBL, since it depends 
on the unknown redshift of the source. The spectral slope is very 
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Fig. 6. Zoom on the RXTE/PCA, XRT and UVOT light curves 
around the flare observed with XRT. Horizontal error bars in- 
dicate the duration of the time intervals over which fluxes were 
averaged. A flux decrease can be seen in the UVOT UV data that 
accompanies the decline in the X-ray flux after the flare. 



steep, which might be at least partially due to absorption on the 
EBL. 

In the high-energy y-ray band, the rather flat ferm/-LAT 
spectrum, extracted over a period of more than four months of 
relatively constant flux (cf. Section I3.1l l. suggests a wide peak, 
followed by the steep spectrum in the VHE range. 

In the X-ray range, the combined spectrum including Swift 
XRT and RXTE PCA data for MJD 55184 is included in the 
SED. These data represent an average flux level between the 
observed peak and the lowest observed flux in the Swift XRT 
light curve. Swift UVOT data and optical data from the ATOM 
telescope are also included for this night. The optical and X-ray 
data have been corrected for Galactic absorption as described 
in Section 14.11 Archival radio d a ta fro m sever al catalogues 
(PKSCAT90, Wright & Otrupcekl (Il990l): PMN, [Wright et al. 



1994); SUMSS, Mauchetal. (2003); CRATES, Healey et al. 
(2007); AT20G. lMurphv et al.1 (120 1 ol) ; ATPMN. fMcConnel et ail 
d2012in have been added to the SED. 

The SED clearly shows the double-bumped structure that is 
characteristic of BL Lac objects. The peak frequency of the high- 
energy bump is relatively well constrained with the information 
from the ,Ferm;-LAT and H.E.S.S. data, if a certain redshift of 
the source and EBL model are assumed. The detection of a flare 



in X-rays and a correlated variation in the UV flux suggests a 
common origin for the emission in these wavelength bands (at 
least above the R band). Given the high flux level in the optical 
and UV range — as compared to the lower X-ray flux — and 
the rather hard X-ray spectrum, the synchrotron peak must be 
located not far above or within the UV range. It is possible that 
at the relatively long wavelengths seen in the R band, another 
component dominates the optical flux. There might be a con- 
tribution from the host galaxy, but given the featureless optical 
spectrum and the non-detection of an extension of the source so 
far, a non-thermal component seems more likely. An extended 
jet might be responsible for the emission at low energies, while 
the emission at higher energies could stem from a more compact 
region inside the jet. This would explain that variability is seen 
in the X-ray and UV range, but not (significantly) in the opti- 
cal bands. The existence of different radiative components in the 
optical and high-energy b ands in blazars is frequently found (see 
Abramowski et al. (2012) for a recent example). 

4.3. Upper limit on the redshift of PKS 0447-439 

The combined ferm/-LAT and H.E.S.S. data set is used to derive 
an upper limit on the redshift of the source. No assumptions are 
made on the source emission characteristics, other than exclud- 
ing spectral upturns beyond the fen«/-LAT energy band, which 
are not expected in general and would be theoretically difficult 
to account for. 

One case of a spectral upturn at high energi es might 
have been seen during a flare of Mrk 501 dNeronovl 
ISemikoz & Tavlorl) . iLefa. Aharonian & Rieger I (1201 ll) interpret 
this hypothetical spectral upturn as resulting from the combi- 
nation of several emission regions in a leptonic scenario with 
Maxwell-like particle disbributions. The observational evidence 
for this feature, however, is weak. For PKS 0447-439 , there is 
no indication of an upturn in the ferm/-LAT data. 

Spectral upturns might also occur in a scenario where the 
y-ray flux in the ferm/-LAT band suffers from internal absorp- 
tion in the source, e.g. on e mission from the Broad Line Region 
(BLR) dSenturk et al.l201 UlPoutanen & Stern |2012|) or from the 
accretion disk or torus (e.g. iDone a & Protheroe 200 3J), and re - 
covers in the VHE band (see also IZacharopoulou et *aT1l201 ll) . 
However, the absence of line features in the optical spectra from 
BL Lac objects in general suggests that emission from the BLR 
is negligible in these sources. In the case of PKS 0447-439, the 
available optical spectra, as discussed in Section Q] do not show 
any significant features, except for absorption lines from the 
Earth's atmosphere. The y-ray data from PKS 0447-439 are best 
described by power laws, with no indication for absorption on 
internal photon fields. 

The fen«/-LAT power-law spectrum, constructed over the 
integration period between MJD 55050 to MJD 55185, is extrap- 
olated to the highest energies to serve as an upper limit on the 
unabsorb ed flux in the VHE band. Following th e method pro- 
posed by iGeor ganonoulos. Fi nke & Revesl (120101) . the ratio be- 
tween the VHE flux and the extrapolated Ferm/-LAT flux pro- 
vides an upper limit on the optical depth, as a function of en- 
ergy, due to EBL attenuation. A comparison of this upper limit 
with the predictions of a given EBL model as a function of red- 
shift can be translated into an upper limit on the redshift of 
the source. This method has been used to constrain the density 
of the EBL with data from Fermi-LAT a nd MAGIC from the 
FSRQ PKS 1222+21 (lAleksic et al.ll201 lb . A similar approach 
was applied to ferm/-LAT and VHE data from the BL Lac ob- 
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Fig. 7. Spectral energy distribution of PKS 0447-439 showing archival radio data from NED and SIMBAD/VizieR (black squares; 
details are given in the text), optical and UV data from ATOM (blue filled circles) and Swift UVOT (red filled circles), X-ray data 
from Swift XRT and RXTE PCA (magenta squares), y-ray data from Fermi-LP^ (black triangles and grey bow-tie) and VHE y-ray 
data from H.E.S.S. (blue filled circles and light blue bow-tie). The models shown here correspond to a redshift of 0.2 (black solid 
line; deabsorbed spectrum shown by thin black dashed line), 0.3 (red dashed line; deabsorbed spectrum shown by thin red dashed 
line) and 0.4 (blue dash-dotte d line; deabsorbed spect rum shown by thin blue dashed line). For the graph on the left panel, the EBL 
model uses the description bv lFranceschini et"aD(l2008|). On the right panel, only the high-energy part of the SED is shown, with an 
EBL model following the description by Kneisk e & Dole I (l2010h . which presents a very low EBL level. The inlays show a zoom on 
the X-ray range and VHE range of the SED. 



ject PG 1553+1 13 to d erive an upper limit on its unknown red- 
shift dAbdoet al.ll20Tol) . 



Intrinsic softening of the spectrum beyond the ferm/-LAT 
band, which might well occur in PKS 0447-439, would only 
lead to a lower value for the inferred redshift and would thus 
not be in contradiction with the derived upper limit. A prob- 
lem in the interpretation of the upper limit would only arise 
in non-standard emission scenarios, in which absorption on 
the EBL is partially avoided. This is the case for scenar- 
ios (e.g. lEssev et al.ll201 it iMurase et al.ll2012[: [A haronia n et al] 
I2012t iTavlor et al.ll201 ll lAharonian et al.ll2002l) , where the ob- 
served y-rays are interpreted as secondary photons produced 
outside of the source in interactions of primary photons or 
protons with the cosmological microwave background and the 
EBL. Another possibility is given by scenarios involving ex- 
otic particles or interactions, such as ph oton-axion oscillations 
dDe Angelis, Galanti & Roncadellill2011l) . 



the upper limit on the optical depth r max is given by 



Modifying Equation [Tl from lAleksic et"aT1 (1201 ll) to include 
a correction term for the systematic error on the H.E.S.S. flux, 



i" ma x(£) = In 



F- mt (E) 



(l-a)-(F ohs (E)-l.64-AF obs (E)) 



(1) 



Here Fi nt (E) is the intrinsic flux, assumed to be a simple ex- 
trapolation of the flux measured with ferm/-LAT, F \, S (E) is the 
VHE flux measured with H.E.S.S., with a statistical uncertainty 
of AF b s (E), and (1 - a) is the factor accounting for the sys- 
tematic error on F^iE). The systematic error is a measure of 
the uncertainty in the flux and is expressed as a fraction of the 
measured flux. For the analysis cuts used here and for the mea- 
sured spectral index, the factor a equals about 0.26, as deter- 
mined from detailed Monte Carlo simulations of the instrumen- 
tal performance, corresponding to a ~ 26% systematic error on 
the flux averaged between 0.3 TeV and 2 TeV This is taken into 
account in a conservative way, in the sense that a is the maxi- 
mum factor by which the flux could be overestimated. 

As shown in Fig. [8] an upper limit of z<0.45 can be put on 
the redshift of the source at the 95% confidence leve l (CL), when 
using the EBL model by iFranceschini etakl (120081) . which gen- 
erally provides a good description of the observed VHE spectra. 
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z=0.58 (KD 10) 
=0.59 (low, G 09) 
z=0.45 (F 08) 
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Energy [TeV] 




0.5 0.55 
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Fig. 8. Left panel: upper limit at 95 % CL on the optical depth as a function of energy (solid line) as derived from Equation Q] 
The optical depth distribution derived with different EBL models at specific redshifts is shown for comparison (dash-dotted line: 
Kneiske and Dole 2010; dashed line: Gilmore et al. 2009, "low" model; dotted line: Franceschini et al. 2008. Right panel: the 
minimum difference 5 m i„,E between the maximum optical depth derived here and the optical depths predicted by the same EBL 
models is shown as a function of the redshift. Eo is the energy for which, at a given redshift, the difference between the derived 
and predicted optical depths is minimized. The upper limit on the redshift, for each EBL model under consideration, corresponds to 



- , min,Eo 



=0. 



Applying an estimation of the lo wer limit of the EBL level, as 
provided by the "low" model bv iGilmore et al.l (120091) and by 
the model by Kneiske & Dole] ( 1201 Oh . a more conservative up- 
per limit of z<0.59 and z<0.58 is derived, respectively. 

A detailed discussion of the inclusion of statistical and sys- 
tematic uncertainties of the FermZ-LAT spectrum in the upper 
limit, not included in Equation [1] is given in Appendix [Al It is 
shown that this approximation is indeed justified in the present 
case. 

An estimate of the redshift of PKS 0447-439 
of z= 0.20±0.05 has rec en tly been suggested by 



Prandini, Bonn oli & Tavecchio d201ll). based on prelimi- 
nary data from H.E.S.S. (|ZecheL^|20111) and on an empirical 



relation bet ween observed photo n indices in the Fermi-LAT and 
VHE band (Prand ini et alJl2010l) . Such an estimate is based on 
the assumption that all detected VHE blazars share very similar 
intrinsic spectral properties. The most conservative upper limit 
of z<0.59 derived in the present work, on the other hand, pro- 
vides a firm limit on the redshift of the source, which takes into 
account statistical and systematic uncertainties in the spectral 
measurements and depends only on minimal assumptions on the 
spectral shape and the validity of current EBL models. 

5. Interpretation with an SSC Model 

The interpretation of the multiwavelength spectrum of 
PKS 0447-439 is complicated by the fact that the redshift of the 
source is unknown, but it is possible to study the variation of the 
model parameters for different redshift values within the derived 
lower and upperlimits. The modelling presented here, based on a 
simple SSC model, aims at investigating two issues: the physical 
parameters characterizing the emission region and the maximum 
allowable redshift of the source that can be accommodated under 
the assumption of a simple SSC scenario. 

The stationary homogeneous one-z one SSC model used for 
this s tudy is based on the code by Kat arzvnski. Sol & Kus I 
(2001). Synchrotron and SSC emission stem from a single popu- 
lation of relativistic electrons and positrons in a spherical plasma 



blob inside the jet. This blob is characterized by its radius R, 
Doppler factor 5 and homogeneous, tangled magnetic field B. 
The energy distribution of the relativistic particles is assumed to 
follow a broken power law with indices n\, rt2, minimum and 
maximum Lorentz factors y m [ n and y max , a break at a Lorentz 
factor ybi and a normalisation K\ at a Lorentz factor of y — 1 . 

When applying the usual light travel time arguments, an up- 
per limit can be derived on the ratio of the size of the emission 
region (blob) R and its relativistic bulk Doppler factor 6 from the 
observed variabililty time scale St b S : 



R6~ 



1 +z 



■ Af, 



obs ■ 



(2) 



The flux variability time scale in the X-ray data of about one 
day constrains this ratio to R6 l < 2.2xl0 15 cm for the lower 
redshift limit of z=0.176 (cf. SectionQ]) and to R6' 1 < 1.6xl0 15 
cm for the upper limit of z=0.59, derived in Section 1431 These 
values are not unusual for BL Lac objects. 

The synchrotron and SSC peaks are not directly detected, 
but the peak frequencies can be constrained from the available 
data. The high flux level and hard spectrum in the optical and 
UV band, compared to the relatively low flux in the X-ray band, 
exclude a synchrotron peak position that is too far above the UV 
data and imply a peak frequency between 10 15 and 10 16 Hz in 
the observer's frame. 

In the Thomson limit, the Lorentz factor at the break in the 
electron spectrum (yy) is defined by the ratio of the observed 
synchrotron and SSC peak frequencies, v s and v c , respectively 
dTavecchio et al.lll998l) : 



Tbr 



3v_ 

4v« 



1/2 



(3) 



This implies a value of y^ between 10 4 and 10 5 , when assuming 
an SSC peak position between 10 24 and 10 26 Hz. 

A minimum Doppler factor is give n by the opacity condi- 
tion according to lTavecchio et al.l (fl998) and lDondi & Ghisellinil 
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d!995h : 



6> 



— d 2 (l 



F(v ) 



Afobs 



1/(4+2/?) 



(4) 



Here di is the luminosity distance of the source, F(vq) the flux 
at the frequency of the target photons, /3 the spectral index of 
the target photons and cr T the Thomson cross-section. For the 
minimum redshift of z=0.176, this leads to a minimum Doppler 
factor of S >14. Assuming a redshift of z=0.3, shown in the 
following to be still acceptable for a one-zone SSC description, 
this value would increase to 6 >18. 

In the first two scenarios to be investigated, the redshift is 
fixed at the lower limit of z=0.176 and at z=0.20, respectively. 
Good representations of the SED can be found for the param- 
eters listed in Table [5] (first two columns); the modelled SED 
corresponding to z=0.20 (black solid line) is shown in Fig. |7] 
for ab sorption on the EBL wi th the model by Frances chini et al.l 
( 2008) and with the model bv lKneiske & Dole U2010h . Emission 
at the lowest frequencies in the radio band is ascribed to the ex- 
tended radio lobes (cf. Section [TJ and is not modelled here. 

To arrive at these solutions, a low magnetic field of B =0.02 G 
and relatively high bulk Doppler factors of <5=35 and £=37.5 are 
assumed. The low magnetic field is necessary to account for the 
low synchrotron peak frequency. Increasing B, while keeping all 
other parameters fixed, shifts the synchrotron peak to higher fre- 
quencies and prevents a good representation of the X-ray slope. 
This can be compensated by lowering the Doppler factor accord- 
ingly, but then the overall flux level in the model decreases. To 
boost the flux level, the electron density K has to be increased, 
leading to a dominance of the SSC peak that cannot be compen- 
sated with an increase of the radius of the emission region, due to 
the constraint from Equation|2] Satisfactory solutions with lower 
Doppler factors can be found down to 6 ~30. For smaller 6, it 
becomes increasingly difficult to reproduce the similar flux level 
in the synchrotron and SSC bumps, while staying in agreement 
with the other constraints from the data. 

As is often seen in SSC scenarios for BL Lac objects, the en- 
ergy density in electrons and positrons and the magnetic energy 
density are clearly out of equipartition, with u e /UB ~ 70. 

The minimum Lorentz factor of the electron distribution y m ; n 
has to be of the order of a few 100 due to the upper flux limit that 
is provided by the radio data, although it has to be pointed out 
that these data are not simultaneous. The break energy of the 
electron distribution depends on the synchrotron and SSC peak 
frequencies as discussed above. 

An interpretation of the SED has also been attempted for 
higher redshifts, up to z=0.4 (cf. Table and Fig. |7). Only the 
bulk Doppler factor 5 has been increased to compensate for the 
decreasing flux and spectral shift with higher redshift. With in- 
creasing z, absorption on the EBL becomes more severe and the 
model has a tendency to cut off below the H.E.S.S. spectrum. 
Trying to correct for this effect with an increasing break energy 
in the electron spectrum does not improve the result sufficiently 
and is limited by the position of the synchrotron peak. The model 
does not provide a satisfactory representation of the H.E.S.S. 
data for redshifts larger than about 0.35 for the EBL model by 
Francesc hini et al.l d2008l) and f or redshifts larger t han ab out 0.4 
for the minimum EBL model bv lKneiske & Dole I (|2010), as can 
be seen from Fig Q In addition, the Doppler factor required to 
compensate for the effect of increasing redshift, becomes very 
large for the solutions at high z. 

It can also be seen that, in the one-zone SSC scenario, the 
SSC peak frequency has to be located between 10 24 and 10 25 
Hz. 
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400 






Ybt 






2.5 x 10 4 






/max 






9.0 x 10 5 






u e [erg crrr 3 ] 






1.6 x 10" 3 






ii B [erg cm -3 ] 






2.3 x 10" 5 







Table 5. Parameters used for stationary SSC modelling of the 
SED of PKS 0447-439 for different assumptions on the redshift 
of the source. A description of the parameters is given in the text. 



6. Summary 

The distant HBL PKS 0447-439 has been detected in the VHE 
band for the first time, with a significance of 15.2cr for 13.5 
hours of live time. Emission at these energies follows a steep 
spectrum with a photon index of 3.89+0.37 and shows marginal 
evidence of variability on a monthly scale. During the H.E.S.S. 
observations, rapid flaring was detected in the X-ray range, ac- 
companied by a flux variation in the UV range, at a time scale 
of about one day. This variation is not reflected in the Fermi- 
LAT and VHE energy bands, but might well be hidden in the 
statistical uncertainties and by gaps in the coverage. This puts a 
limit on the size of the emission region, which is comparable to 
observations from other HBL. 

A comparison of the Fermi-LAT spectrum with the H.E.S.S. 
spectrum leads to a conservative upper limit on the redshift of 
th e source (z<0 .59 at 9 5% CL for the lower limit EBL model 
bv lGilmore et al.l (12009)), taking into account statistical and sys- 
tematic uncertainties. This is a firm upper limit, under the as- 
sumptions that the intrinsic flux does not exhibit any upturns 
above the Fermi-LAT band and that current models provide a 
correct estimate of the minimum density of the EBL. 

Assuming a simple SSC scenario to explain the SED of 
PKS 0447-439 from the optical to the VHE range, the physi- 
cal parameters of the emission region and underlying particle 
distribution were investigated for different redshifts. Good solu- 
tions at low redshifts require a low magnetic field and large size 
of the emission region, close to the upper limit from variability 
constraints. With a one-zone, homogeneous SSC model, a satis- 
factory description of the SED from the optical to the VHE range 
does not seem feasible for redshifts above about 0.4. 
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Appendix A: Treatment of the uncertainties in the 
Fermi-LAT flux for a derivation of the redshift 
upper limit 

Apart from the approximative treatment of the systematic error 
in the H.E.S.S. flux, EquationQ]provides an exact expression for 
the 95% C.L. upper limit on the maximum optical depth r mflA , if 
the intrinsic spectrum of the source is known precisely. However, 
if the statistical uncertainty on the Fermi-LAI spectrum is to 
be taken into account, an evaluation of the confidence interval 
is complicated by the fact that T max is a function of the ratio 
of two Gaussian variables, the observed fluxes in the HE and 
VHE range, and is thus itself not a Gaussian variable. Moreover, 
the ratio of two Gaussian variables does not possess moments to 
any order greater or equal to one. Therefore, confidence intervals 
must rely on percentiles. 

As will be shown in the following, the confidence interval 
[D- n ,D +n ] at the n-cr confidence level for the ratio of two inde- 
pendent Gaussian variables with mean a, /3 and standard devia- 
tion cr a , crp is given by the expression: 



a ± 



D + 



(a> 2 +P' 2 -n 2 r- 



(a' 1 +l3* 2 -n 2 )1 



na-p 



(A.l) 



Here, a* = a/cr a , (3* = p/crp and n is the n-cr parameter 
defining the confidence level y (see Table lA. ll for a few examples 
of y as a function of n). This can be seen as follows: suppose one 
is interested in the ratio D — y/x of two independent Gaussian 
random variables x and y : x ~ N(J3, cri), y ~ N(a, cr 2 a ). The 
joint probability distribution function of the couple (x, y) is 



f(x,y) = — exp{--[( 

Z7T 2 LV (Tp 



1 \(X- yS\2 ,y - 



(A.2) 



Next, the change of variables x = pcrp cos and y = pcr a sin 
is performed. Knowing 0, one knows the ratio D, therefore one 
seeks the probability distribution function of 0, which is 



i 

2n~ 



pe 



dp . 



with 



1 tan 0o = 



The integral in ( IA.3I > can be evaluated easily. It yields 

1 r 1 P- 

yfhiR cos(0 - o )e"^ 2 ^-^(R cos(0 - O ))] , 



(A.3) 



(A.4) 



M<P) 



(A.5) 



where <t> is the cumulative distribution function of a reduced 
Gaussian variable, that is 



O(f) 



•J -co 



1 



V2^ 



du . 



(A.6) 



From now on, two approximations are introduced. First, the 
signal-to-noise ratio R is supposed to be sufficiently high, such 
that e~~- R ~ x 0. Second, the variations of (p about 0o are supposed 
not to be too extreme, such that <f>(R cos(0 - 0o)) as 1 . This last 
condition is certainly fulfilled if all quantities of interest stay 
positive. Under this condition, the trigonometric functions are 



n 


1 1.64 


2 


2.58 


3 


7 


68.3% 90.0% 


95.4% 


99.0% 


99.7% 



Table A.l. Confidence level y as a function of the n-cr parameter 
for typical values of n and y. 



one-to-one relations of their arguments. Equation (IA.5) reduces 
to 

M<f>) = -Lfl cos(0 - ft) expf-i/? 2 sin 2 (0 - ft,)] . (A.7) 

V27T 

If one sets t — R sin(0 - 0o)> me probability distribution function 
(p.d.f.) of t equals 



1 

V2lr" 



ft(t)= -i=eM~\t 2 ), 



(A.8) 



which is the p.d.f. of a reduced Gaussian variable. Since there are 
one to one relations between t, <p and the ratio D, the percentiles 
of D are uniquely determined by the percentiles of t . Next, one 
defines t n , f_„ = -t„ and the confidence level y by 



I 



1 , i 

— exp(- k 2 ) dt = <D(n) - O(-n) = y(n) . (A.9) 
t-„ V2^ 



That is 

t ±n = ±n . 
From this one derives successively 

R sin(0 ± „ - O ) = ±n , 
+n 

4>± n - 0o = arcsin — , 

K 

n 



arctan 



(D, 



<p ±n = 0o ± arcsin - , 
R 

, — = arctan ± arcsin — , 



n 
R 



era ,a eras n 

D ± „ — = tan(arctan ) ± arcsin — ) 

cr a y Pcr a > R 



(A. 10) 

(A. 11) 
(A. 12) 

(A.13) 

(A. 14) 

(A.15) 



Now, with the use of the trigonometric relations: tan(a + b) — 
(tana + tanb)/(l +tanatan/?) and tan arcsin x = xj Vl - x 2 , one 
gets 

P °~l> {R 2 -n 2 )l 



D + „ = 



1 + I- 



' (R 2 -n 2 )1 



(A. 16) 



-/r 



Introducing a* = a/cr a and j3* = filcrp, such that R 2 = . 
one finds Equation dA.U . 

The effect of including the uncertainties on the Fermi-LAT 
spectrum on the resulting distribution of T max (E) can be seen 
from Fig. IA.1I where the exact result derived here, i.e. D + „ 
with n=1.64, is compared to the distribution arrived at when ne- 
glecting the statistical uncertainty in the Fermi-LAT spectrum, 
and with an approach where H.E.S.S. and Fermi-LAT errors are 
treated in the linear approximation, i.e. by simple error prop- 
agation. Systematic uncertainties have been neglected here. It 
should be noted that in the energy range of interest for the deriva- 
tion of an upper redshift limit, around 1 TeV, all three approaches 
provide very similar results. 

For the curves shown in Fig. IA.2I an estimate of the sys- 
tematic uncertainty in the H.E.S.S. and Fermi-LAT fluxes has 
been added as a bias, as described in Sec. I4.3l for H.E.S.S. For 
Fermi-LAT, a systematic bias of 10% has been assumed. In the 
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Fig. A.l. Distribution of T max vs. energy for three different treat- 
ments of the statistical errors (systematic errors are neglected), 
red triangles: only statistical errors in the H.E.S.S. spectrum are 
taken into account; black circles: statistical errors in the Fermi- 
LAT and H.E.S.S. fluxes are included in the linear approxima- 
tion; blue open squares: statistical errors from Fermi-LKI and 
H.E.S.S. are included using the exact analytic solution derived 
here. 



LAT flux and the inclusion of the statistical uncertainties in the 
Fermz'-LAT flux have a negligible effect. 
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Fig. A.2. Distribution of T max vs. energy including systematic er- 
rors as flux biases, green circles: identical to the curve in Fig. [8] 
following Eq. [1] red triangles: same as above, but with a sys- 
tematic error of 10% added to the Ferm/-LAT spectrum (cf. 
Sec. I3.lt : blue open circles: statistical errors from Ferm/-LAT 
and H.E.S.S. are included using the exact solution derived here, 
systematic errors have been added to both fluxes as well. 



energy range of interest, the exact solution is in good agreement 
with the curve used to derive the upper limit on the redshift of 
the source (cf. Fig. [8j. The systematic bias added to the Fermi- 
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